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Context. The Abell 222 and 223 clusters are located at an average redshift z~0.21 and are separated by 0.26 deg. Signatures of mergers have 
been previously found in these clusters, both in X-rays and at optical wavelengths, thus motivating our study. In X-rays, they are relatively 
bright, and Abell 223 shows a double structure. A filament has also been detected between the clusters both at optical and X-ray wavelengths. 
Aims. We analyse the optical properties of these two clusters based on deep imaging in two bands, derive their galaxy luminosity functions 
(GLFs) and correlate these properties with X-ray characteristics derived from XMM-Newton data. 

Methods. The optical part of our study is based on archive images obtained with the CFHT Megaprime/Megacam camera, covering a total 
region of about 1 deg 2 , or 12.3x12.3 Mpc 2 at a redshift of 0.21. The X-ray analysis is based on archive XMM-Newton images. 
Results. The GLFs of Abell 222 in the g' and r' bands are well fit by a Schechter function; the GLF is steeper in r' than in g'. For Abell 223, the 
GLFs in both bands require a second component at bright magnitudes, added to a Schechter function; they are similar in both bands. The Serna 
& Gerbal method allows to separate well the two clusters. No obvious filamentary structures are detected at very large scales around the clusters, 
but a third cluster at the same redshift, Abell 209, is located at a projected distance of 19.2 Mpc. X-ray temperature and metallicity maps reveal 
that the temperature and metallicity of the X-ray gas are quite homogeneous in Abell 222, while they are very perturbed in Abell 223. 
Conclusions. The Abell 222/ Abell 223 system is complex. The two clusters that form this structure present very different dynamical states. 
Abell 222 is a smaller, less massive and almost isothermal cluster. On the other hand, Abell 223 is more massive and has most probably been 
crossed by a subcluster on its way to the northeast. As a consequence, the temperature distribution is very inhomogeneous. Signs of recent 
interactions are also detected in the optical data where this cluster shows a "perturbed" GLF. In summary, the multiwavelength analyses of 
Abell 222 and Abell 223 are used to investigate the connection between the ICM and the cluster galaxy properties in an interacting system. 
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1 . Introduction faint-end slopes is a good way to trace the history of the faint 

galaxy population and the influence of merging events on this 

Galaxy evolution is known to be influenced by environmental population . As summarize d for example in Table A. 1 of Boue 

effects, which are particularly strong in galaxy clusters, all the et al (2008)j M& slope can strongly vary from one cluster t0 

more when they are undergoing merging events. The analysis anotheri and from one region t0 another in a given duster 
of galaxy luminosity functions (GLF), and in particular of their 

Superseding the dynamical analyses of Proust et al. (2000) 

* Based on observations obtained with MegaPrime/MegaCam, a and Girardi & Mezzetti (2001), Dietrich et al. (2002) have ob- 

joint project of CFHT and CEA/DAPNIA, at the Canada-France- tained a redshift catalogue for 183 galaxies in the region of 

Hawaii Telescope (CFHT) which is operated by the National Research the cluster pair Abell 222/223, out of which 153 belong to one 

Council (NRC) of Canada, the Institut National des Sciences de of the two c i usters . T^jr ana l ys i s was a l so based on photom- 

r Univers of the Centre National de la Recherche Scientifique (CNRS) etry in ^ y an( j R bands They have estimated reds hifts of 

of France, and the University of Hawaii. This work has also made use 02126±00008 and 0.2079+0.0008, and velocity dispersions 

of data products produced at TERAPIX and the Canadian Astronomy „ , _, . , , ___ , . „ . , ,, „„„ , . , ,, „„„ 

_ „ c , _, , _ .. _ . T of 1014 and 1032 km s , for Abell 222 and Abell 223 re- 
Data Centre as part ot the Canada-France-Hawaii Telescope Legacy 

Survey, a collaborative project of NRC and CNRS. The X-ray analysis specttvely. As indicated by Dietrich et al. (2002) these velocity 

is based on XMM-Newton archive data. This research has made use of dispersions are somewhat higher than those expected from the 

the NASA/IPAC Extragalactic Database (NED) which is operated by X-ray luminosities. Abell 222 appears relatively relaxed, with 

the Jet Propulsion Laboratory, California Institute of Technology, un- a galaxy velocity histogram consistent with a gaussian. On the 

der contract with the National Aeronautics and Space Administration, other hand, the adaptive kernel galaxy density map of Abell 223 
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shows two peaks, though surprisingly the Dressler-Schectman 
test (Dressier & Schectman 1988) does not find these two peaks 
and the DIP statistics cannot reject the null hypothesis of a uni- 
modal distribution, implying that a Gaussian parent population 
cannot be excluded. Dietrich et al. (2002) derived respective 
mass to light ratios of 202+43 and 149+33 h 70 M Q /L Q and lu- 
minosity functions in the R band down to absolute magnitudes 
of -19.5, and found that a bridge of galaxies seems to exist be- 
tween the two clusters. Note that such bridges of galaxies join- 
ing two clusters are still quite rare, another example at a compa- 
rable redshift being the filament detected between Abell 1763 
and Abell 1770 by Fadda et al. (2008). 

A weak lensing analysis of the Abell 222/223 system was 
performed by Dietrich et al. (2005), who found evidence for a 
possible dark matter filament connecting both clusters. They 
interpret the difference in redshift between Abell 222 and 
Abell 223 as implying a separation along the line of sight of 
15+3 h^ 1 Mpc. Their weak lensing surface mass density con- 
tours again show evidence for a double structure in Abell 223. 

The existence of a 1.2 Mpc wide filament between both 
clusters was also found in X-rays by Werner et al. (2008). 

We present here a muti-wavelength optical and X-ray anal- 
ysis of the clusters Abell 222 and Abell 223. The optical part 
is based on imaging in the g' and r' bands obtained with the 
CFHT Megaprime/Megacam camera and retrieved from the 
Megapipe image stacking pipeline at CADC (Gwyn 2009). 
As seen in Table Q] the r' band data are about 0.55 magni- 
tude deeper than the g' band data (by taking an average color 
g' - r' — 0.85), and will be given priority in the morphological 
analysis; the g' data will be used to select bright cluster mem- 
bers from color-magnitude diagrams, and also because the g' 
band is more sensitive to star formation than the r' band. These 
images are deeper than those of Dietrich by at least 0.5 mag 
(our raw galaxy counts start decreasing for r' > 24.5 while 
Dietrich et al. say that their counts in the R band stop following 
a power law for R> 24) , and therefore allow to probe the GLF 
to somewhat fainter magnitudes. Our X-ray analysis is based 
on archive XMM-Newton data. 

We will consider hereafter that the two clusters have a com- 
mon redshift of 0.21, for which Wight's cosmology calculator 
(Wright 2006) gives a distance of 1035 Mpc, a spatial scale of 
3.427 kpc/arcsec and a distance modulus of 40.07 (assuming a 
flat ACDM cosmology with H = 70 km s" 1 Mpc" 1 , Q M = 0.3 
and Q. A = 0.7). 

The paper is organised as follows. We describe our opti- 
cal analysis in Section 2, and results concerning the galaxy lu- 
minosity function in Section 3. A search for substructures in 
both clusters based on the Serna-Gerbal method is described in 
Section 4. The X-ray data analysis and results, including tem- 
perature and metallicity maps, are presented in Section 5. An 
overall picture of this cluster pair is drawn in Section 6. 

2. Optical data and analysis 

2.1. The optical data 

We have retrieved from the CADC Megapipe archive (Gwyn 
2009) the reduced and stacked images in the g' and r' 



Table 1. Summary of the observations. 



Filter 



Number of coadded images 
Exposure time (s) 
Seeing (arcsec) 
Limiting magnitude (5cr) 



5 8 

2000 4000 

0.81 0.73 

26.7 26.4 



bands (namely G002.024.493-12.783.G and G002.024.493- 
12. 783. R) and give a few details on the observations in 
Table Q] Observations were made at the CFHT with the 
Megaprime/Megacam camera, which has a pixel size of 0.1 86 x 
0.186 arcsec 2 . 

We did not use the catalogues available for these images, 
because they were made without masking the surroundings of 
bright stars, so we preferred to build masks first, then to extract 
sources with SExtractor (Bertin & Arnouts 1996). Due to the 
dithering pattern, the total area covered by the r' image was 
slightly larger than the Megacam 1 deg 2 field: 1.196 deg 2 . 

Objects were detected and measured in the full r' im- 
age, then measured in the g' image in double image mode. 
Magnitudes are in the AB system. The objects located in the 
masked regions were then taken out of the catalogue, leading 
to a final catalogue of 223,414 objects. All these objects have 
measured r' magnitudes, and 214,669 also have measured g' 
magnitudes. 

Since the r' catalogue is deeper, and the seeing is also bet- 
ter in this band (see Table [TJ, we will perform our star-galaxy 
separation in the r' band. 

2.2. Star-galaxy separation 

In order to separate stars from galaxies, we plotted the maxi- 
mum surface brightness fi max in the r' band as a function of r' . 
The result is shown in Fig.Q] 

The best fit to the star sequence visible on Fig.[T]calculated 
for 17.25 < r' < 22 is /x max = 0.996 * r' - 0.247, with standard 
deviations on the slope and constant of 0.002 and 0.047 re- 
spectively. The point-source (hereafter called "star") sequence 
is clearly visible for r' < 22, with the star saturation show- 
ing well for r' < 17.25. We will define galaxies as the ob- 
jects: with Hmax(r') > 17.25 for r' < 17.25, and as the objects 
above the line of equation fimax = 0.996 * r' — 0.247 + 0.3, 
or Umax - 0.996 * r' + 0.053. Stars will be defined as all the 
other objects (see Fig. Q]). The small cloud of points observed 
in Fig. Q] under the star sequence is in fact defects, but repre- 
sents less than 2% of the number of stars. We thus obtain a star 
and a galaxy catalogue. 

As a check to see up to what magnitude we could trust our 
star-galaxy separation, we retrieved the star catalogue from the 
Besangon model for our Galaxy (Robin et al. 2003) in a 1 deg 2 
region centered on the position of the image analysed here. 
Such a catalogue is in AB magnitudes (as ours) and is corrected 
for extinction. In order for it to be directly comparable to our 
star catalogue, we corrected our star catalogue (and our galaxy 
catalogue as well, for later purposes) for extinction: 0.071 mag 
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Objects detected in r 





Fig. 1. Central surface brightness in the r' band as a function 
of r' magnitude. The horizontal and oblique full red lines sepa- 
rate the galaxies (above these lines) from the stars (below these 
lines). The oblique dotted red line shows the limit for stars com- 
puted for 17.25 < r' < 22 (see text). The two vertical lines cor- 
respond to r' - 17.25 where the stars stop being saturated, and 
r' = 22, where the fit to calculate the star sequence was limited 
(see text). 



in g' and 0.052 mag in r' (as derived from the Schlegel et al. 
1998 maps). 

The r' magnitude histogram of the objects classified as stars 
in our r' image roughly agrees with the Besanijon star catalogue 
for r' <20, though in some bins our counts tend to be some- 
what higher than those of the Besan^on model even for bright 
magnitudes. For r' >20, we detect more stars than predicted by 
the Besanc,on model (the difference is about 35% at r' ~ 21). 
This could imply that our star-galaxy separation is not valid for 
r' >20; however, as suggested by Fig.[T] this is a rather bright 
limit, as confirmed by simulations adapted to match compara- 
ble data, where we found that the star-galaxy separation de- 
termined with this method was reliable at least up to r' - 21 
(Boue et al. 2008). In reality, the disagreement between our 
star counts and those predicted by the Besancon model is most 
probably due to the fact that our field is located in the direction 
of one of the densest regions of the Sagittarius stream (Ibata et 
al. 2001, Yanny et al. 2009). It is therefore normal to observe 
more stars than predicted by the Besangon model, which does 
not take this stream into account. Note that there may also be a 
small contribution due to quasars and AGN, since their number 
at r' ~ 20.5 is expected to be about 20 per square degree (see 
Richards et al. 2006, figure 13). 

For r' > 20, we will conservatively consider that the star- 
galaxy separation may be wrong, and we will compute galaxy 
counts by counting the total number of objects (galaxies plus 
stars) per bin of 0.5 mag, and considering that the number 
of galaxies is equal to the total number of objects minus the 
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Fig. 2. Positions of the objects (stars and galaxies) in the 
Abell 222 (magenta) and Abell 223 (cyan) catalogues. The ob- 
jects used to estimate background counts (see text) are shown in 
blue for the north rectangle and in red for the annulus surround- 
ing the clusters. Note that the figure covers 1.2 x 1.2 deg 2 , and 
is slightly larger than the images, which cover 1 .04 x 1 . 1 5 deg 2 . 



number of stars predicted in each bin by the Besangon model. 
As mentioned below, we have checked that the galaxy number 
counts in the 20 < r' < 22 range estimated by both meth- 
ods (i.e. our star-galaxy separation and the subtraction of the 
Besangon counts to the total number of objects) are consistent 
within error bars. 

Since the clusters are quite distant, they do not cover the 
whole field, so we extracted from the star and galaxy catalogues 
two catalogues as large as possible corresponding to the two 
clusters. We defined the J2000.0 positions of Abell 222 and 
Abell 223 as coinciding with the brightest galaxy of each clus- 
ter, i.e. 24.3921,-12.9912 and 24.5099,-12.7575 respectively 
(in degrees). Note that these are not exactly the cluster posi- 
tions given by the NED database. The maximum possible ra- 
dius to obtain independent catalogues for the two clusters was 
0.1290 deg, or 1.6 Mpc at a redshift of 0.21. 

For the Abell 222 and Abell 223 clusters separately, we 
thus obtained three complementary catalogues (with g' and r' 
magnitudes): objects classified as galaxies (classification valid 
at least for r' < 20), objects classified as stars, and a complete 
catalogue of galaxies+stars which will be used for r' > 20. The 
positions of the galaxies in the regions of the two clusters are 
shown in Fig. [2] 



2.3. Catalogue completeness 

The completeness of the catalogue is estimated by simulations. 
For this, we add "artificial stars" (i.e. 2-dimensional Gaussian 
profiles with the same Full- Width at Half Maximum as the av- 
erage image Point Spread Function) of different magnitudes to 
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Fig. 3. Point source completeness as a function of magnitude in 
percentages for Abell 222 (top) and Abell 223 (bottom) in g' 
(left) and r' (right) for point-like objects (see text). 

the CCD images and then attempt to recover them by running 
SExtractor again with the same parameters used for object de- 
tection and classification on the original images. In this way, 
the completeness is measured on the original images. 

In practice, we extract from the full field of view two 
subimages, each 1500 x 1500 pixels 2 , corresponding to the po- 
sitions of the Abell 222 and Abell 223 clusters on the image. 

In each subfield, and for each 0.5 magnitude bin between 
r' - 20 and 27, we generate and add to the image one star 
that we then try to detect with SExtractor, assuming the same 
parameters as previously. This process is repeated 100 times 
for each of the two fields and bands. 

Such simulations give a completeness percentage for stars. 
This is obviously an upper limit for the completeness level for 
galaxies, since stars are easier to detect than galaxies. However, 
we have shown in a previous paper that this method gives a 
good estimate of the completeness for normal galaxies if we 
apply a shift of ~ 0.5 mag (see Adami et al. 2006). Results are 
shown in Fig. [3] 

From these simulations, and taking into account the fact 
that results are worse by ~ 0.5 mag for mean galaxy popula- 
tions than for stars, we can consider that our galaxy catalogue 
is complete to better than 80% for g' < 25 and r' < 24.5. 

2.4. Galaxy counts 

The surfaces covered by the Abell 222 and 223 catalogues (af- 
ter excluding masked regions) are 0.0429 deg 2 and 0.0475 deg 2 
respectively. Galaxy counts were computed in bins of 0.5 mag 
normalized to a surface of 1 deg 2 . 

For r' < 20, galaxy counts were derived directly by com- 
puting histograms of the numbers of galaxies in the Abell 222 
and Abell 223 catalogues. For r' > 20, we built for each cluster 
histograms of the total numbers of objects (galaxies+stars) and 
obtained galaxy counts by subtracting the numbers of stars pre- 
dicted by the Besancon model. The resulting galaxy counts will 



be used in the next section to derive the GLFs for both clusters 
in both bands. 

As a test, we considered the galaxy counts in the four 0.5 
magnitude bins for 20 < r' < 22 computed for both clusters 
with the two methods (i.e. first method: considering that the 
star-galaxy separation is valid, and second method: considering 
the total number of objects (galaxies+stars) and subtracting the 
number of stars predicted by the Besangon model to obtain the 
number of galaxies). In all cases, the differences are smaller 
than 15% (and in most cases they are only a few percent). 
Therefore, the limit (between r' - 20 and 22) at which we 
consider that our star-galaxy separation is not valid any more 
will not strongly influence our results. 

Note that no k-correction was applied to the galaxy magni- 
tudes. 

3. Results: colour-magnitude diagrams and galaxy 
luminosity functions 

In order to compute the galaxy luminosity functions of the two 
clusters, it is necessary to subtract to the total galaxy counts 
the number counts corresponding to the contamination by the 
foreground and background galaxies. 

For galaxies brighter than r' - 20 we will select galax- 
ies with a high probability to belong to the clusters by draw- 
ing colour-magnitude diagrams and selecting galaxies located 
close to this relation. A few spirals may be missed in this way, 
but their number in any case is expected to be small, as ex- 
plained at the end of Sect. 3.1 (also see e.g. Adami et al. 1998). 
For galaxies fainter than r' - 20 we will subtract galaxy counts 
statistically. 

3.1. Colour-magnitude diagrams 

The mean values for g' - r' are 0.85 and 0.87 for Abell 222 
and Abell 223. The g' - r' vs. r' colour-magnitude diagrams 
are shown in Fig. |4]for the two clusters. A sequence is well 
defined for galaxies in the magnitude range 18 < r' < 22 in 
both clusters. We computed the best fit to the g' - r' vs. r' re- 
lations in this magnitude range by applying a linear regression 
(separately for the two clusters, since their colour-magnitude 
relations may not be absolutely identical). We then eliminated 
the galaxies located more than 3<x away from this relation and 
computed the g' - r' vs. r' relations again. 

The equations of the colour-magnitude relations are found 
to be: g'-r' = -0.0545r'+2.091 with 3cr = 0.21 for Abell 222, 
and g'-r' = -0.0474/ + 1.945 with 3cr = 0.21 for Abell 223. 

For r' < 20, we will consider hereafter that all the galaxies 
located within 3cr of these relations (i.e. between the two black 
lines of Fig. 3]i belong to the clusters. With this condition, there 
are 141 and 144 galaxies belonging to Abell 222 and Abell 223 
respectively with r 1 < 20. 

We also plot in Fig.|4]the galaxies belonging to the clusters 
according to their spectroscopic redshifts, i.e. galaxies in the 
[0.195, 0.215] redshift interval (see Sect. 4). We can see that 
their positions agree well with the colour-magnitude selection. 

For both clusters, we have estimated the number of blue 
cluster galaxies lost by selecting galaxies in the red sequence 
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Fig. 4. (g' — r') vs. r' colour-magnitude diagrams for Abell 222 
(top) and Abell 223 (bottom) for objects classified as galax- 
ies from the \i max -magnitude relation. The vertical dashed lines 
indicate the magnitude interval where the colour-magnitude re- 
lation was computed. The long oblique line shows the colour- 
magnitude relation, while the two short oblique lines indicate 
an interval of +3<x around the colour-magnitude relation. The 
filled circles show the galaxies belonging to the clusters ac- 
cording to their spectroscopic redshifts. 



interval in the following way. First, we computed histograms 
of numbers of galaxies in the red sequence and below the red 
sequence in bins of 1 absolute magnitude in the r' band. The 
bins of interest here are M r < = -21.5 and -20.5 (see Fig. 0), 
roughly corresponding to r' - 18.5 and 19.5. Then we esti- 
mated the number of foreground galaxies expected. Since the 
comoving volume at z=0.21 is 2.622 Gpc 3 (Wright 2006), and 
each of our clusters covers an area of 0.0522 deg 2 on the sky, 
the volume in the direction of each cluster is 3318 Mpc 3 . By 
using the R band luminosity function by Ilbert et al. (2005) in 
the [0.05-0.20] redshift bin (see their Fig. 6 and Table 1), we 
find that the percentages of "lost" galaxies are of the order of 
10%-25% forM r - = -20.5 and of 10% for M r , = -21.5. 
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Fig. 5. Galaxy counts in the g' (left) and r' (right) bands for 
magnitudes r' > 20 where the background must be subtracted 
statistically, in logarithmic scale. The counts in Abell 222 and 
223 are drawn in magenta and cyan respectively. The black 
dashed lines show the galaxy counts from the CFHTLS Deep 
survey, and the full green line shows the average of the Deep 1, 
Deep 2 and Deep 3 field counts. The blue and red lines cor- 
respond to the "local" backgrounds extracted from the same 
images as our clusters (red: in the annulus between the two cir- 
cles of Fig. [2] blue: in the north rectangle, see text and Fig. [2j. 
Error bars are Poissonian and are not plotted for clarity. 

3.2. Comparison field 

In order to perform a statistical subtraction of the background 
contribution for r' > 20, we considered several possibilities. 

The most obvious solution would be to extract the back- 
ground in an annulus surrounding the clusters. The largest outer 
circle that can be extracted, centered on the middle position be- 
tween the two clusters (24.4510, -12.8744 J2000.0 in degrees), 
has a radius of 0.4472 deg, or 5.55 Mpc. We can then take an 
inner circle of radius 0.4028 deg, or 5 Mpc. This annulus is 
drawn in Fig. [2] However, the number counts in this circle are 
higher than the CFHLTS Deep survey counts (see below), prob- 
ably because the annulus is too close to the clusters and still 
includes cluster galaxies, so we will not use these counts as a 
background. 

The Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) has been taken with the same telescope, 
camera and filters as the data we are analyzing. The 
Deep survey explores a solid angle of 4 x 1 deg 2 
of the deep Universe, in four independent fields 
(http://www.cfht.hawaii.edu/Science/CFHLS/). Observations 
are carried out in five filters {u*,g',r',i' and z') providing 
catalogs of sources that are 80% complete up to ('^=26.0 (see 
http://terapix.iap.fr/cplt/oldSite/Descart/CFHTLS-T0005-Releas"eTpd^ l. 
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We only consider the Deep survey here, because we are in- 
terested in the GLFs down to faint magnitudes. For the four 
deep fields, we computed the galaxy number counts in bins 
of 0.5 mag, in g' and r', normalized to a surface of 1 deg 2 . 
The galaxy counts in these four fields are drawn in Fig. [5] 
As discussed before (see e.g. Boue et al. 2008 and references 
therein), the counts in these four fields somewhat differ, due 
to cosmic variance, and in particular the counts in the Deep 4 
field are higher than in the other three. We will therefore take 
as background galaxy counts the average of the Deep 1, Deep 2 
and Deep 3 field counts. 

For comparison, we also extracted galaxy counts from a 
region assumed to be representative of the background in our 
image. This region was chosen to be a rectangle north of the 
clusters, in an area devoid of bright stars. It is shown as a blue 
rectangle in Fig.|2]and covers an effective area of 0.06825 deg 2 . 
The corresponding counts (hereafter the "local" background 
counts) appear to be notably lower than the CFHTLS-Deep 
counts (Fig. [5]); this may be the case if this region corresponds 
to a void located between the large scale filaments that con- 
verge towards the clusters. 

Although we consider that the CFHTLS-Deep counts prob- 
ably represent better the galaxy background counts, we will 
also fit the GLFs with these "local" background counts, keep- 
ing in mind the fact that these "local" counts are likely to over- 
estimate the GLFs thus derived for Abell 222 and 223. This 
will illustrate the difficulty to determine GLF parameters un- 
ambiguously. 

We can note from Fig. [5] that galaxy counts are notably 
higher in Abell 223 than in Abell 222, implying that the for- 
mer cluster is richer than the latter. 

3.3. Galaxy luminosity functions 

The Galaxy Luminosity Functions (GLFs) of Abell 222 and 
Abell 223 were calculated in bins of 0.5 mag and normalized 
to 1 deg 2 , as described above. We subtracted the background 
contribution using as background galaxy counts: 1) the average 
of the three Deep 1, Deep 2 and Deep 3 fields; 2) the "local" 
counts in the north blue rectangle of Fig. [2] 

The GLFs are displayed for Abell 222 and Abell 223 in 
Figs.|6]and|7]respectively. The error bars drawn in these figures 
were taken to be 4 times the Poissonian errors on galaxy counts, 
as suggested by detailed simulations previously performed by 
our team for similar data (see Boue et al. 2008, Fig. 5). 

The GLFs (as a function of absolute magnitude) were fit by 
a Schechter function: 



S(M) = 0.4 In 10 c 



'y a+1 e- 



withy= io a4 < M *- M >. 

The parameters of the Schechter function fits of the GLFs 
are respectively given in Tables [2] and [3] for the two different 
background subtractions: the average of the CFHTLS Deep 
fields and the "local" background. The absolute magnitude 
range considered is indicated for each fit. The fits are drawn 
in Figs. [6] and [7] for Abell 222 and Abell 223 respectively (with 
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Fig. 6. Galaxy luminosity functions for Abell 222 in the g' (top) 
and r' (bottom) bands, in logarithmic scale. The pooints corre- 
spond to the subtraction of the background counts taken from 
the CFHTLS. Error bars are 4 times the Poissonian errors on 
galaxy counts (see text). The best Schechter function fits are 
drawn in red for the subtraction of background counts taken 
from the CFHTLS and in blue for "local" background counts 
(see text). 



background galaxy counts taken from the CFHTLS Deep field 
counts). 

If we compare the results given in Tables [2] and [3] we see 
that the cD* and M* parameters depend little on the background 
galaxy counts chosen. On the other hand, some variations are 
found in the faint end slope a, as expected since it is at faint 
magnitudes that the background galaxy subtraction has a strong 
influence. 

For Abell 222, the GLF is quite well fit by a Schechter func- 
tion in the r' band, with a faint end slope a = -1.43 (the same 
for both background subtractions). On the other hand, the GLF 
is surprisingly flat in g' (slope a = -0.60 or -0.94, depending 
on background subtraction). 

A Schechter function is obviously not sufficient to fit the 
GLFs of Abell 223, both in g' and r'\ Abell 223 has more 
very bright galaxies than Abell 222, and a second component at 
bright magnitudes is obviously required. At faint magnitudes, 
the GLF slopes in g' and r' are comparable for Abell 223. They 
vary by small quantities (0.12 and 0.05 respectively) between 
one galaxy background subtraction and the other. However, 
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Table 2. Schechter parameters for galaxy luminosity functions (background galaxy counts taken from the CFHTLS Deep field 
counts). 



Cluster 


Filter 


Range 


O* 


M* 


a 


Abell 222 


ft 


[-22.5,-16.0] 


2891 ± 134 


-20.04 ± 0.06 


-0.60 ± 0.05 




r' 


[-23.5,-16.0] 


937 ± 64 


-22.01 ± 0.07 


-1.43 ±0.01 


Abell 223 


ft 


[-25.5,-15.5] 


1017 ±72 


-21.05 ±0.06 


-1.30 ±0.02 




r' 


[-26.0,-15.0] 


1107 ±62 


-22.05 ± 0.06 


-1.29 ±0.01 



Table 3. Schechter parameters for galaxy luminosity functions ("local" background galaxy counts). 



Cluster 


Filter 


Range 


O* 


M* 


a 


Abell 222 


ft 


[-22.0,-16.0] 


2269 ± 133 


-20.32 ± 0.06 


-0.94 ± 0.03 




r' 


[-23.0,-16.0] 


1132 ±63 


-21.85 ±0.06 


-1.42 ±0.01 


Abell 223 


ft 


[-25.0,-15.0] 


819 ±55 


-21.20 ±0.05 


-1.42 ±0.01 




r' 


[-26.0,-15.0] 


1042 ± 58 


-22.10 ±0.06 


-1.34 ±0.01 
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Fig. 7. Same as Fig.|6]for Abell 223. 



these slopes are only indicative, since as seen in Fig.[7]the GLFs 
show strong wiggles for absolute magnitudes between -20 and 
-17 in r' and between -20 and -16 in g' . Since Abell 223 is 
itself a double cluster (Dietrich et al. 2002), it is not surprising 
to find that it has a "perturbed" GLF This is also the case for 
example for the Coma cluster, known to have a two-component 
GLF (Biviano et al. 1995). 

If we compare the GLFs of the two clusters, we can see 
that they are more or less comparable in the r' band (except 



for the wiggles in Abell 223), while in the g' band the GLF is 
notably steeper in Abell 223 than in Abell 222. An explanation 
for this can be that both clusters have comparable old galaxy 
populations, but that the number of faint star forming galaxies 
is higher in Abell 223, where star formation can be triggered 
by recent interactions. 

If we compare these GLFs to those found by other au- 
thors, we can note that the bright part of the GLF Schechter 
fits (y < 18, or M,- < -18.9) for both clusters are very similar 
in shape to the GLFs recently obtained e.g. by Andreon et al. 
(2008); these authors analyzed the GLFs of a sample of clusters 
at various redshifts, limited to relatively bright absolute magni- 
tudes: M v < -19. 

The a slopes of the faint ends of the GLFs derived here 
(except for that of Abell 222 in the g' band) are within the 
broad range of values estimated by previous authors for dif- 
ferent clusters, cluster regions and photometric bands (see e.g. 
the compilation in Table Al of Boue et al. 2008). 

Note that Abell 222 and Abell 223 are at redshift 0.21, and 
few GLFs are available for clusters at such redshifts. Andreon 
et al. (2005) found somewhat shallower slopes for three clusters 
at redshifts ~ 0.3, but in the K band, so the comparison with our 
results is not straightforward. 

3.4. Morphological segregation 

We applied a new tool developed in SExtractor that calculates 
for each galaxy the respective fluxes in the bulge (spheroid) 
and disk. This new experimental SExtractor feature fits to 
each galaxy a two-dimensional model comprised of a de 
Vaucouleurs spheroid (bulge) and an exponential disk. Briefly, 
the fitting process is very similar to that of the GalFit pack- 
age (Peng et al. 2002) and is based on a modified Levenberg- 
Marquardt minimisation algorithm. The model is convolved 
with a supersampled model of the local Point Spread Function 
(PSF), and downsampled to the final image resolution. The PSF 
model used in the fits was derived with the PSFEx software 
(Bertin et al. 2010, in preparation) from a selection of 5436 
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Fig. 8. Fraction of disk-dominated galaxies as a function of ra- Fig. 10. Average r' magnitude as a function of radius for 
dius for Abell 222 (in magenta) and Abell 223 (in cyan). Abell 222 (in magenta) and Abell 223 (in cyan). 




n/N is: 



Fig. 9. Fraction of disk-dominated galaxies in sectors of 45 deg 
within each cluster, with the same colours as in the previous 
figure. Angular sectors are numbered from 1 to 8 clockwise 
from east (see Fig. IA. U . 



(g r ) I 6922 (r') point source images. The PSF variations were 
fit using a 6th degree polynomial of x and y image coordinates. 

The model fitting was carried out independently in the g' 
and r' bands. The bulge to total ratios discussed here were ex- 
tracted from the red channel, limited to r' - 22, since beyond 
this magnitude results may become unreliable. 

We applied this tool to look for evidence for morphologi- 
cal segregation, by computing the fraction of galaxies with a 
spheroid to total flux smaller than 0.3, as a function of radius 
(i.e. in concentric circles of step 0.02 deg). The result is plotted 
in Fig. [8] In each bin, lo - errors were computed as follows: if n 
is the number of galaxies having a spheroid to total flux smaller 
than 0.3 and N the total number of galaxies, we assume that n 
and (N-n) are independent variables; in this case, the error on 



(n/N) y/(l/n-l/N). 



We can see that there is a general increase of the fraction 
of disk-dominated galaxies with radius as expected (see e.g. 
Biviano et al. 1997). 

In order to see if there was a morphological segregation that 
could be linked to the presence of cosmological filaments, we 
computed that same fraction in angles of 45 deg all around each 
cluster. The result is shown in Fig. [9] where angular sectors are 
numbered from 1 to 8 clockwise from east (see Fig. IA. lb . No 
strong trend is found in this plot within error bars, as expected 
if there is no obvious cosmological filament linking the clusters 
with the surrounding cosmic web. We can note however that for 
both clusters there is a dip in the fraction of bulge-dominated 
galaxies in sector 6, which corresponds to the line joining the 
two clusters; the filament already detected by Werner et al. 
(2008) using wavelet-decomposition in the [0.5-2] keV energy 
band is located along this same direction. In this zone the gas 
is not yet hot enough to influence the galaxies. If we look at 
figure 9 from Lagana et al. (2009), 2 keV is the lower limit 
where ram-pressure plays an important role in galaxy clusters. 
Thus we would expect to detect a trend only if the filament was 
hotter than 2 keV. 

Another possibility for this dip in the fraction of bulge- 
dominated galaxies is a higher number of blue galaxies in 
which star formation was triggered due to ongoing collisions. 
Poggianti et al. (2004) argued that there is a striking corre- 
lation between the positions of the young and strong post- 
starburst galaxies and substructure in the hot intra-cluster 
medium (ICM) identified from XMM-Newton data, with these 
galaxies lying close to the edges of two infalling substructures. 
This result suggests that the interaction with the dense ICM 
could be responsible for the triggering of the star formation. 
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Fig. 11. Positions of the galaxies with measured redshifts and 
magnitudes in the regions covered by Abell 222 and Abell 223 
(black crosses). The galaxies belonging to the two main grav- 
itationally bound systems identified with the SG method are 
shown in colours (see text). 

3.5. Luminosity segregation 

We searched for luminosity segregation in both clusters by se- 
lecting galaxies brighter than r' - 24 and within the red se- 
quence defined in Fig. [4] to increase the probability that they 
belong to the clusters. The mean r' magnitude was calculated in 
each of the radial bins previously defined and results are shown 
in Fig[10] We observe an increase of the mean magnitude with 
radius: it is about 21.3 in the central bin, and then increases 
with radius from values around 22 to 22.6. 

4. Cluster-scale to large-scale structure 

We searched the NED database for galaxies with redshifts 
available in a very large region of 3 deg radius around 
Abell 222/223 and found 980 galaxies. This radius corresponds 
to 37 Mpc at the distance of the studied clusters, i.e. to the typ- 
ical diameter of a "bubble" or of a void in the universe (see 
e.g. Hoyle & Vogeley 2004). This catalogue of galaxies with 
measured redshifts will be used first to search for substructures 
in the zones covered by our clusters, and second to search for 
structures (such as filaments) at a larger scale. 

4.1. Search for gravitationally bound structures 

Within the regions of our images covered by Abell 222 and 
Abell 223, 246 galaxies have measured redshifts, and 210 
galaxies have both measured redshifts and magnitudes. Out of 
these 210 galaxies, 173 have redshifts in the [0.18,0.24] inter- 
val. The positions of the 210 galaxies with measured redshifts 
and magnitudes are displayed in Fig. QT| It is interesting to 
note that, as traced by these galaxies, both clusters appear to 



Fig. 12. Redshift histograms for all the galaxies with measured 
redshifts in the regions covered by Abell 222 and Abell 223 (in 
black), and for the galaxies in the gravitationally bound struc- 
tures of Abell 222 and Abell 223, in magenta and cyan respec- 
tively (the positions of these galaxies are shown in Fig.fTTTi with 
the same colour coding. 



be strongly elongated, with elongation position angles that dif- 
fer from one cluster to the other. This strongly suggests that the 
clusters are not in dynamical equilibrium, in agreement with 
previous studies (see above, and also Zabludoff et al. 1995, 
Roettiger et al. 1996, Boschin et al. 2004, Girardi et al. 2006). 

The corresponding redshift histogram, zoomed on the clus- 
ter redshift interval, is shown in Fig. [12] 

We applied to the catalogue of 210 galaxies with redshifts 
and magnitudes the Serna & Gerbal (1996) algorithm (here- 
after SG), which allows to search for substructures, separate 
galaxies forming gravitationally bound structures, and esti- 
mate the masses of these substructures. The dendogram thus 
obtained is shown in Fig. [13] It clearly shows two dynami- 
cally distinct substructures, corresponding to the two clusters: 
Abell 222 in the top half and Abell 223 in the bottom half. 
The two main gravitationally bound structures corresponding 
to these two clusters respectively include 55 and 64 group 
members. Their corresponding velocity dispersions are 1014 
and 1170 km s~ : respectively. The virial masses derived for 
these two substructures (assuming M/L=400 solar units) are 
1.9 10 11 and 4.2 10 13 M . 

For Abell 222, and probably also for Abell 223, 
these masses are obviously underestimated, particularly for 
Abell 222, since such a value does not even reach that of a 
cD galaxy. We have tried to analyze the two clusters sepa- 
rately, but this modifies the resulting masses only slightly. The 
fact that we find too small cluster masses can be due to sev- 
eral reasons. First, the redshift sample we are using is taken 
from the NED data base, and is obviously not complete, even 
at bright magnitudes. We have estimated the completeness of 
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Fig. 13. Dendogram showing how Abell 222 and Abell 223 can 
be dynamically separated when applying the Serna & Gerbal 
algorithm. The abscissa is the negative binding energy while 
the catalogue numbers of the various galaxies are shown along 
the ordinate (the top of the dendogram is truncated, but would 
only show a constant gravitational binding energy). 



the redshift catalogue by comparing the numbers of galaxies 
with redshifts in the [0.18,0.24] redshift interval to the total 
number of galaxies, in magnitude bins of 0.5 in the r' band. We 
find that for 17 < r' < 19.5 the completeness is around 50%, 
then rapidly decreases for r' > 19.5, confirming that the cluster 
masses based on such a redshift catalogue are probably unreli- 
able. Second, as pointed out by the referee, the SG method re- 
lies on the unlikely assumption that galaxies in clusters are the 
mass carriers, and additionally requires an assumption about 
the M/L ratio, here taken to be M/L=400; this value may be 
reasonable for a cluster taken as a whole, but is too large for in- 
dividual galaxies (the dark matter halos of individual galaxies 
would be overlapping in the dense central regions of clusters, 
and strong mass segregation would result by dynamical fric- 
tion, which is not observed). 

Better estimates of the cluster masses can be obtained by 
applying other methods. For example, the relation between 
mass and velocity dispersion computed by Biviano et al. (2006) 
leads to masses of 1.2 10 15 and 1.4 10 15 M (with H = 70 
km s -1 Mpc~') for Abell 222 and Abell 223 respectively, based 
on respective velocity dispersions of 1014 and 1070 km s . 
Note however that since the clusters are not virialized these val- 
ues are probably overestimates. For example, Takizawa et al. 
(2010) have recently shown from numerical simulations that 



for merging systems where the subcluster has a mass larger 
than one fourth that of the large cluster, the virial mass can be 
larger than twice the real mass. This most probably applies at 
least to the Abell 223 cluster. 

We can also apply the mass-temperature relation estimated 
from X-ray data by Arnaud et al. (2005), with the overall tem- 
peratures given in Table 5; this relation gives masses M200 of 
3.8 10 14 and 4.7 10 14 M for Abell 222 and Abell 223 respec- 
tively (also with Ho = 70 km s _1 Mpc~'). These values may be 
closer to the real values than those estimated from the galaxy 
velocity dispersions, since e.g. Takizawa et al. (2010) found 
that for merging clusters X-ray derived masses were usually 
more reliable than virial mass estimations. However, as ex- 
plained for example by Nagai et al. (2007), although the to- 
tal ICM mass can be measured quite accurately (to better than 
~6%) in clusters, the hydrostatic estimate of the gravitationally 
bound mass is biased low by about 5-20% throughout the virial 
region, primarily due to additional pressure support provided 
by subsonic bulk motions in the ICM, ubiquitous in our simu- 
lations even in relaxed systems. Another source of error is the 
fact that the cluster temperatures, particularly that of Abell 223 
(see Fig. [T8l are not homogeneous. For example, Rasia et al. 
(2004) have shown that an incomplete thermalisation of the gas 
due to an incomplete virialisation of the cluster could lead to an 
underestimate of the mass. On the other hand, if the X-ray gas 
is heated by a collision, the mass could be overestimated. 

We can also note that Mamon (1993) has calculated the bias 
on the M/L ratio when the virial theorem is wrongly applied to 
a group which is still evolving dynamically (see his figure 1 1). 
He found that in some cases the mass could be underestimated 
by a factor reaching 100. Therefore the masses estimated for 
our clusters assuming that they are virialized are clearly not 
valid. 

Unfortunately, the SG method does not allow to separate 
the two substructures in Abell 223 discovered by Dietrich et al. 
(2002) based on galaxy density contours (see their figure 6). 

Note that Abell 223 is probably close to the intersection 
of two or more filaments, since it is massive and substruc- 
tured, and its GLF is comparable of that e.g. of Coma (see 
Section 3.3). On the other hand, Abell 222 is probably only 
following its path along a filament. 



4.2. Large-scale structure 

Out of the 980 galaxies found in NED in a 3 deg radius re- 
gion, 310 have redshifts in the [0.18,0.24] interval (i.e. about 
±9000 km s -1 around the mean cluster velocities). Their spa- 
tial distribution is shown in Fig. [14] A cluster is visible towards 
the southwest: Abell 209. It has a redshift similar to that of 
Abell 222 and Abell 223 and is located about 19.2 Mpc away 
from Abell 222 in projection on the sky. No significant filamen- 
tary structure is detected towards the Abell 222/223 clusters. 
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Fig. 14. Distribution of the galaxies with available redshifts in 
a large region surrounding the clusters. Blue points: all galax- 
ies with measured redshifts, red points: galaxies with redshifts 
in the [0.18,0.24] interval. The two clusters under study are cir- 
cled. A third cluster appears to the southwest: Abell 209. The 
physical size of this map is 44.7 x 44.7 Mpc 2 . 



The background was normalized using a spectrum obtained in 
an annulus (between 12.5-14 arcmin) where the cluster emis- 
sion is no longer detected. In one of the pointings, the A223 
cluster falls exactly inside the annulus where we determine the 
background contribution. We thus masked this region, avoid- 
ing any cluster contamination. In Fig. [16] we show a compari- 
son between the observed and Read backgrounds for the three 
detectors, as well as the residuals between both determinations. 




Fig. 15. Annulus where the background was extracted for sub- 
traction. 



5. X-ray analysis 

5.1. Data reduction 

The pair of clusters Abell 222/223 was observed with XMM- 
Newton in two different pointings (revolutions 1378 and 1380) 
with a total exposure time of 144 ks. We used data from all 
EPIC cameras (MOS1, MOS2 and pn). The data were reduced 
with the XMM-Newton Science Analysis System (SAS) v8.0 
and calibration database with all updates available prior to 
November 2009. The initial data screening was applied us- 
ing recommended sets of event patterns, 0-12 and 0-4 for the 
MOS and PN cameras, respectively. 

The light curves are not constant and large variations in in- 
tensity are visible (flares). To improve the signal-to-noise ratio 
we discarded periods of flares and the cleaned light curves in 
the energy range of [1-10] keV exhibited stable mean count 
rates; exposure times are given in Table [4] We considered 
events inside the field of view (FOV) and excluded all bad pix- 
els. 

Table 4. Information on X-ray data 



Name 



Filter 



^exp (kS) 

MOS1 MOS2 



n H (10 2U cm- 2 ) 



pn 



Abell 222 
Abell 223 



Thinl 
Thinl 



11.814 

23.479 



18.005 
26.069 



10.226 
10.612 



2.26 
2.26 



Table 5. Overall X-ray temperatures and metallicities 



Cluster 



''200 

(Mpc) 



kT 
(keV) 



Z 
(solar) 



Abell 222 
Abell 223 



1.28: 
1.55: 



0.11 
0.15 



3.77: 
4.38: 



0.15 
0.16 



0.23: 
0.23: 



0.06 
0.08 



The background was taken into account by extracting 
MOS1, MOS2 and pn spectra from the publicly available EPIC 
blank sky templates of Andy Read (Read & Ponman 2003). 



In order to determine the global properties of these clusters, 
we adopted the virial radii r2oo determined by the weak lens- 
ing analysis (Dietrich et al. 2005). These values, together with 
the mean temperatures and metallicities are given in Table [5] 
They were calculated fixing the hydrogen column density to its 
galactic value of «h = 1-56 x 10 20 cm -2 . 



5.2. Spectrally measured 2D X-ray maps 

We performed quantitative studies using X-ray spectrally mea- 
sured 2D maps to derive global properties of these clusters. 
These maps were made in a grid; for each spatial bin we set a 
minimum count number of 900 (after background subtraction). 
For the spectral fits, we used XSPEC version 11.0.1 (Arnaud 
1996) and modeled the obtained spectra with a MEKAL sin- 
gle temperature plasma emission model (bremsstrahlung + line 
emission, Kaastra & Mewe 1993, Liedahl et al. 1995). The free 
parameters are the X-ray temperature (kT) and the metal abun- 
dance (metallicity). Spectral fits were made in the energy inter- 
val of [0.7-8.0] keV with the hydrogen column density fixed 
at the galactic value (see Table @), estimated with the nH task 
of FTOOLS (based on Dickey & Lockman 1990). 
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Fig. 17. Merged event maps from MOS1 (left panel), MOS2 (central panel) and pn (right panel) showing the regions where the 
2D maps were obtained. 
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spectra: 2 MOS1, 2 MOS2 and 2 PN). The best fit value is then 
attributed to the central pixel. 
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Fig. 16. Comparison of the observed and Read backgrounds for 
MOS1 (top panel), MOS2 (middle panel) and pn (low panel). 
The black points correspond to the observed data, the red points 
to Read's scaled background and the green points to the resid- 
uals. 



We compute the effective area files (ARFs) and the re- 
sponse matrices (RMFs) for each region in the grid. This pro- 
cedure (already presented in Durret et al. 2005 and Lagana et 
al. 2008) allows us to perform a reliable spectral analysis in 
each spatial bin, in order to derive high precision metallicity 
and temperature maps (see Fig. [T~ ). since we simultaneously 
fit all three instruments for both pointings (that is, we coadd 6 
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Fig. 18. X-ray temperature map for Abell 222 (bottom) and 
Abell 223 (top). The color-bar indicates the temperature in keV. 
The corresponding error map is shown in Fig. IB. II 

The resulting X-ray temperature and metallicity maps are 
displayed in Figs. [18] and [19] Maps of the errors on these two 
maps are given in the Appendix (Figs. lB~T1 and lB.2t . 

The study of the thermal structure of the intra-cluster 
medium (ICM) provides a very interesting record of the dy- 
namical processes that clusters of galaxies have experienced 
during their formation and evolution. The temperature (and 
pressure) distribution of the ICM gives us important insight 
into the process of galaxy cluster merging and on the dis- 
sipation of the merger energy in form of turbulent motion. 
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Fig. 19. X-ray metallicity map for Abell 222 (bottom) and 
Abell 223 (top). The color-bar indicates the metallicity in solar 
units. The corresponding error map is shown in Fig. 



Metallicity maps can be regarded as a record of the integral 
yield of all the different stars that have released their metals 
through supernova explosions or winds during cluster evolu- 
tion. 

From the temperature maps we can see that neither cluster 
presents a defined cool-core as a sign of relaxation. However, 
Abell 222 is fairly isothermal when compared to Abell 223. 
The overall temperature structure in Abell 223 shows the 
clear imprint of the recent interaction(s) this cluster has been 
through. The cluster exhibits a remarkable blob of hotter gas at 
its southwest extremity, a colder region to the southeast, and 
a hotter zone crossing the cluster along a roughly east-west 
direction. The metallicity map of Abell 222 is quite homoge- 
neous, while that of Abell 223 is very inhomogeneous, with a 
metallicity enhancement in a zone which does not completely 
coincide with the colder southeast region and a low metallic- 
ity region extending along a zone curving from the east to the 
west-northwest. 

5.3. Interpretation 

As has been shown by Kapferer et al. (2006, see his figure 
10), the enhancement of the metallicity from Z=O.2Z (in the 
region of the shock) to Z=O.4Z (in the surrounding region) 
takes 3 Gyr (assuming their model B). This means that the gas 
has not yet had time to mix and distribute the metals in the 
region where the shock is still propagating. As can be seen in 
the larger variations in metallicity that Abell 223 presents when 
compared to Abell 222, the metallicity of merging systems is 
not smoothly distributed. It is likely that Abell 223 has recently 



been crossed by a smaller cluster, which now appears at the 
north east tip of the maps. 

From the metallicity map we also observe that the central 
region of Abell 223 presents a metallicity enhancement when 
compared to the surrounding medium, though not highly cen- 
trally peaked. At this point, it is important to highlight that 
Abell 223 has more very bright galaxies than Abell 222. The 
results of a former work (Lagana et al. 2009) suggested that 
BCGs alone cannot produce the observed metallicity excess in 
the central region (De Grandi & Molendi 2001, De Grandi et 
al. 2004). Thus, it is likely that the cool-core has been disrupted 
by the sub-cluster collision but the metals injected into the ICM 
by very bright galaxies of Abell 223 have not moved that much 
away from the central region yet. This can occur if one assumes 
ram-pressure stripping with tidal disruption of these very bright 
galaxies near the center as a mechanism for metal injection. In 
this case, we would expect Abell 223 to have a central metal- 
licity higher than Abell 222, and this is indeed the case. 
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Fig. 20. Surface brightness profiles of the five brightest galaxies 
in Abell 223. The "cD" is in black, the very bright galaxy near 
the cD is in red, and the three other galaxies are somewhat more 
distant from the cluster center. The cD profile decreases notably 
slower with radius that the profiles of the other galaxies. 



It is interesting to note that Abell 223 has in fact two BCGs. 
The galaxy that we consider as being in the cluster center looks 
like a cD, with several smaller satellite galaxies. The surface 
brightness profiles of the brightest galaxies of Abell 223 are 
displayed in Fig. [20] We can see that the surface brightness 
profile of the "cD" decreases slowlier than those of the other 
bright galaxies. It can therefore be considered as a dominant 
galaxy, since it also shows clear signs of an associated popula- 
tion of fainter galaxies in the likely process of being accreted. 
A second nearby elliptical galaxy which is even brighter could 
be the central galaxy of an accreted group. This is reminiscent 
of the Coma cluster where the "cD" (NGC 4874) is not quite as 
bright as the brightest galaxy (NGC 4889), which is believed 
to have been accreted after the cluster was formed. 
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In the case of Abell 222, the situation is even more compli- 
cated, since there are three brightest cluster galaxies, aligned 
along the direction seen in Fig.QT|and having exactly the same 
magnitude, the middle one being surrounded by several small 
satellite galaxies that give it the appearance of a cD. This is an- 
other argument in favour of the idea that Abell 222 is not fully 
relaxed either. 

6. Discussion and Conclusions 

Understanding the interplay between galaxies and ICM is fun- 
damental for further investigations of galaxy clusters. In this 
way, an important step is made with multiwavelength studies 
of these structures. 

In order to understand better the connection between ICM 
and galaxies in merging systems we combined both X-ray 
and optical data to examine in detail the pair of clusters 
Abell 222/ Abell 223. To do so, we combined spectroscopi- 
cally computed 2D-maps with optical data to gain insight on 
the cluster history. We verified that besides their proximity, 
these clusters show very different evolutionary stages that can 
be detected in optical and X-ray analyses. The Abell 223 data 
indicate that this cluster is dynamically perturbed when com- 
pared to its companion. Signs of recent merger(s) are detected 
at both wavelengths, meaning that the gas and the galaxies are 
still out of equilibrium. It is more likely that the signs of merger 
detected in Abell 223 are due to the interaction with the sub- 
cluster. On the other hand, Abell 222 is a smaller cluster closer 
to a stable dynamical state than Abell 223. 

As mentioned in Sect. 15.21 Abell 223 exhibits a remark- 
able blob of hotter gas at its southwest extremity. The gas of 
this region may have been compressed and heated by a recent 
shock that this cluster has been through. Another zone of higher 
temperature is elongated towards the north east-south western 
direction. This is close to the direction that connects the two 
clusters and where a filament of low temperature was already 
detected between these clusters (Werner et al. 2008). This re- 
gion does not appear in our 2D maps because there are not 
enough X-ray counts. However, the X-ray contours overlaid in 
the 2D maps show that the emission goes beyond the limits of 
the maps and towards the region where the filament is located. 

The presence of only a small number of hot blobs and the 
low X-ray temperature of the filament tend to indicate that this 
system is in the first stages of a merger. This hints at a scenario 
in which the two clusters are at the first stages of a collision. 
These clusters are beginning to collide and from the optical 
counterpart analysis, the central region of Abell 223 already 
exhibits signs of this interaction presenting a higher number of 
star forming galaxies. So it seems probable that a small clus- 
ter has just crossed Abell 223 and is heading towards the north 
east. On the other hand, Abell 222 does not present strong evi- 
dence of interaction other than the absence of a cool-core. 

We thus summarize our main results below: 

- From the optical analysis we conclude that Abell 222 is 
smaller and less massive than Abell 223. For Abell 222, 
the GLF is quite well fit by a Schechter function in both 
the r' and g' bands. On the other hand, since Abell 223 is a 



double cluster, it presents a "perturbed" GLF. Abell 223 has 
more bright galaxies than Abell 222 and as a consequence a 
second component is required to fit the GLF at bright mag- 
nitudes. The GLFs of both clusters are comparable in the 
r' band but in the g' band the GLF is notably steeper in 
Abell 223 than in Abell 222. This is explained by a higher 
number of faint star forming galaxies present in Abell 223. 
Star formation is triggered by recent merger events. Such 
mergers disturbed the gas and its imprints can be observed 
in 2D X-ray maps. 

- Both clusters appear to be elongated when traced by their 
cluster galaxy distributions. This is another evidence that 
these clusters are not in equilibrium. However, the X-ray 
filament joining the two clusters does not seem to play a 
very important role in the evolution of cluster galaxies. This 
is probably explained by the low temperature of the gas, 
since ram-pressure is not very efficient below 2 keV. 

- From the spectrally measured 2D X-ray maps we con- 
clude that Abell 222 is fairly isothermal when compared 
to Abell 223. However, neither of these clusters presents 
a cool-core as a sign of relaxation. The overall tempera- 
ture map of Abell 223 is an imprint of the recent dynamical 
interactions this cluster has experienced. The broad varia- 
tion in temperature puts in evidence a large number of sub- 
structures. From the metallicity map, we also observe an 
enhancement in the central region of Abell 223. This can be 
explained by the fact that this cluster has more very bright 
galaxies than Abell 222. Assuming ram-pressure stripping 
with tidal disruption of the bright galaxies near the clus- 
ter center as a mechanism for metal injection, one would 
expect this cluster to have a higher central metallicity than 
Abell 222. The combined analysis of X-ray and optical data 
argues for a scenario in which the two clusters are begin- 
ning to collide. The former cluster presents signs of inter- 
action in both wavelength ranges, and has probably been 
crossed by a small cluster heading north east, while its com- 
panion does not show strong evidence for interaction, ex- 
cept for the absence of a cool-core. 
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Appendix A: Sectors for Fig.|5] 
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Fig. A.l. Positions of the sectors used to draw Fig. [9] 

The positions of the sectors used to draw Fig. [9] are shown 
inFig.lATTI 



Appendix B: Error maps 

The maps of the errors on the X-ray temperature and metal- 
licity maps are displayed in Figs. IB.11 and lB.2l The errors on 
these parameters were directly obtained from the spectral fits 
and they give the deviations in the maps presented in Figs. [T8l 
and[T9] 
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Fig. B.l. Error map on the X-ray temperature. 
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Fig. B.2. Error map on the X-ray metallicity. 



